A class of bridges of iterated integrals of Brownian motion 
related to various boundary value problems involving 
the one-dimensional polyharmonic operator 
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Abstract Let {B{t))t^[Q ^ be the linear Brownian motion and (-''^n(0)te[o,i] be the {n — l)-fold 
integral of Brownian motion, n being a positive integer: 



r (t - sV^i 
Xn{t) - / \ dB{s) for any t E [0, 1] 

Jo (n~iy. 



In this paper we construct several bridges between times and 1 of the process (X„(i))tg[o x] 
involving conditions on the successive derivatives of Xn at times and 1. For this family of bridges, 
we make a correspondance with certain boundary value problems related to the one-dimensional 
polyharmonic operator. We also study the classical problem of prediction. Our results involve 
various Hermite interpolation polynomials. 
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1 Introduction 

Throughout the paper, we shall denote, for any enough differentiable function /, its i-th derivative 
by or d'f/dt\ 

Let {B{t))t^[Q,i] be the linear Brownian motion started at and (/3(i))tG[o,i] be the linear 
Brownian bridge within the time interval [0, 1]: (/3(i))t6[o,i] = {B{t)\B{0) = B{1) = O)tg[o,i]- These 
processes are Gaussian processes with covariance functions 

Cg {s,t) — s At and (s, t) — s At — st. 

For a given continuous function u, the functions and respectively defined on [0, 1] by 

''^si^) — / Cg{s,t)u{s) ds and / c^{s,t)u{s) ds 

Jo Jo 

are the solutions of the respective boundary value problems on [0, 1]: 

and 
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Observe that the differential equations are the same in both cases. Only the boundary condi- 
tions differ. They are Dirichlet-type boundary conditions for Brownian bridge while they are 
Dirichlet/Neumann-type boundary conditions for Brownian motion. 

These well-known connections can be extended to the polyharmonic operator (P"^ /dt^^ where 
n is a positive integer. This latter is associated with the (n — l)-fold integral of Brownian motion 

{Xn{t))te[0,l]- 

Xn{t) = f dB{s) for any t e [0, 1]. 

Jo ("-!)! 

(Notice that all of the derivatives at time naturally vanish: X„(0) = X„_i(0) = • ■ • = ^2(0) = 
Xi{0) — 0.) Indeed, the following facts for instance are known (see, e.g., j9| and (lOj): 

— The covariance fonction of the process (^n(^))te[o.i] coincide with the Green function of the 
boundary value problem 

w(2«) = (_l)"u on [0,1], 
v{0) ^ v'{0) w("-i)(0) = 0, 

«(»)(!) ^ ^...^ v(2«-i)(l) = 0; 

— The covariance fonction of the bridge (X„(t)|X„(l) — 0)fg[o.i] coincide with the Green function 
of the boundary value problem 

w(2») ^ (_l)nif on [0,1], 
v{0) ^v'{0) = ••• = (0) =0, 

w("-i)(l) = t;("+i)(l) = . . . = u(2«-i)(l) = 0; 

— The covariance fonction of the bridge (X„(i)|X„(l) = X„_i(l) = •■. = Xi{l) = 0)t£[o,i] 
coincide with the Green function of the boundary value problem 

w(2«) = (_i)"u on [0,1], 
v{0) = v'{0) = ■■■ = w("-i)(0) = 0, 
t;(l) = v'{l) = ■■■ = w("-i)(l) = 0. 

Observe that the differential equations and the boundary conditions at are the same in all cases. 
Only the boundary conditions at 1 differ. Other boundary value problems can be found in [2] 
and [5]. 

We refer the reader to [3] for a pioneering work dealing with the connections between general 
Gaussian processes and Green functions; see also [1]. We also refer to [2], [2], [8], [H], [12], [15] and 
the references therein for various properties, namely asymptotical study, of the iterated integrals 
of Brownian motion as well as to [3] , [S] , [H] and [H] for interesting applications of these processes 
to statistics. 

The aim of this work is to examine all the possible conditioned processes of {Xn{t))te[o,i\ 
involving different events at time 1: 

(X„(i)|X,,(l) = X,,{1) - • • • = X,„(l) = 0)ie[o,i] 

for a certain number m of events, 1 < m < n, and certain indices ji, j2) ■ ■ ■ ,jm such that 1 < ji < 
j2 < . . . < < and to make the connection with the boundary value problems: 

'w(2«) = (_i)"ii on [0,1], 

< v{0) = v'{0) = ■■■ = w("-i)(0) = 0, (BVP) 
= = . . . w('")(l) = 

for certain indices ii,i2,...,i„ such that < ii < 12 < ••• < in < 2n — 1. Actually, we shall 
see that this connection does not recover all the possible boundary value problems and we shall 
characterize those sets of indices for which such a connection exists. 

The paper is organized as follows. In Section [21 we exhibit the relationships between general 
Gaussian processes and Green functions of certain boundary value problems. In Section [31 we 
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consider the iterated integrals of Brownian motion. In Section we construct several bridges asso- 
ciated with the foregoing processes and depict explicitly their connections with the polyharmonic 
operator together with various boundary conditions. One of the main results is Theorem l4.3l More- 
over, we exhibit several interesting properties of the bridges (Theorems 14.11 and 14. 2 1) and solve the 
prediction problem (Theorems 14.41) . In Section O we illustrate the previous results on the case 
n = 2 related to integrated Brownian motion. Finally, in Section [51 we give a characterization for 
the Green function of the boundary value problem (BVP) to be a covariance function. Another 
one of the main results is Theorem 16.21 



2 Gaussian processes and Green functions 



We consider a n-Markov Gaussian process {X{t))t^[Q ^ evolving on the real line R. By "n-Markov" , 
it is understood that the trajectory 1 1-^ X{t) is n times diffcrcntiable and the n-dimensional process 
{X{t),X'{t), . . . , X("^^^(t))jg[o,i] is a Markov process. Let us introduce the covariance function of 
{X{t))t,£[o,i]- for s,t e [0,1], Cx(s,i) = E[X{s)X{t)]. It is known (see [T]) that the function 
admits the following representation: 

n-l 

c^{s,t) = "^ipk{sAt)Ms'^t) (2.1) 

fc=0 

where (pk, ipk, k G {0, 1, . . . , n — 1}, are certain functions. 

Let VqjVi be linear differential operators of order less than p and let I? be a linear differential 
operator of order p defined by 

i=0 

where ao,ai, .. .,ap are continuous functions on [0, 1]. More precisely, we have for any p times 
differentiable function / and any t G [0, 1], 



(P/)(t) = f]«,(t)/W(t). 

i=0 



Theorem 2.1 Assume that the functions ipk, ipk, k € {0, 1, . . . , n — 1}, are p times differentiable 
and satisfy the following conditions, for a certain constant k: 

r ,W (z) I ] jo tfO<i<P~2, 
J tfi=P~l, 

V^k^V^/jk^O, (2?o¥'fe)(0) = 0, (PiV'fc)(l) =0. (2.3) 
Then, for any continuous function u on [0, 1], the function v defined on [0, 1] by 

v{t) — / c^{s,t)u{s) ds 
Jo 

solves the boundary value problem 

[Vv = a^u, 

Remark 2.1 If the problem ()2.4p is determining, that is if it has a unique solution, then the covari- 
ance function is exactly the Green function of the boundary value problem (|2.4p . 
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Proof In view of (j2.1l) . the function v can be written as 



(i) = ^ tpk{t) (pk{s)u{s) ds + ipk{t) Vfc(s)w(s)ds 



The derivative of v is given by 



k=o I Jt . 



and its second order derivative, since X)fe=o i'Pki/j'k ~ •/'fcV'fc] = 0, by 



fc=0 
n-1 

= E 



fc=0 



'/'fcW / ¥'fc(s)u(s)ds + (p;:(t) / Vfc(s)w(s)ds 



More generally, because of the assumptions (|2.2p . we easily see that, for i G {0, 1, . . . ,p — 1}, 



n— 1 p „ 



/ ^k{sHs)ds + ^''^\t) I MsHs)ds 



and the p-th order derivative of v is given by 



n— 1 r ^ 



E 

fe=0 



^k\t) I MsHs)ds + ^f'{t) I MsHs)d 



Actually, we have proved that, for i G {0,1, . . . ,p — 1} 







(s, t) M(s)ds. 



(2.5) 



Finally, due to 



n— 1 r „ 

Vvit) = E ^^'«(^) / 

fc=o L -^0 



Vipkit) tpkis)u{s) ds + Dtpkit) ■ipkis)u{s) ds 



Kapu{t) — Kapu{t). 



Concerning the boundary value conditions, referring to ()2.3p . we similarly have 

n — 1 ^1 n — 1 „1 

(I)oi^)(0) = E(^o'^fe)(0) / V'fe(s)"(s)ds = 0, (I?i«)(l) = E(^i^fe)(l) / ¥'fe(s)^^(5)ds = 



k=0 



The proof of Theorem 12.11 is finished. □ 



In the two next sections, we construct processes connected to the equation Vv = u subject to 
the boundary value conditions at 0: (X'qw)(0) = for i G {0, 1, . . . , n — 1} and others at 1 that will 
be discussed subsequently, where V and Vq are the differential operators {T> being of order p — 2n) 
defined by 

d^" d* 
^-(-l)"d^' ^"=d^^- 
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3 The (n — l)-fold integral of Brownian motion 



Let (B(t))(g[Q 1] be the linear Brownian motion limited to the time interval [0, 1] and started at 0. 
We introduce the (n — l)-fold integral of Brownian motion: for any t e [0, 1], 



Xnit) = 



(t - s)"-- 
(n-1)! 



■dB{s). 



In particular, Xi — B. The trajectories of (X„(t))jg[o,i] ^ times differentiable and we have 
Xn '^ — Xn^i for < i < rt — 1. Moreover, we have at time the equalities X„(0) = X„_i(0) = 
••• = -^^2(0) — Xi{0) = 0. The process (-'^ri(i))te[o,i] ^ n-Markov Gaussian process since the 
n-dimensional process (X„(t), X„_i(t), . . . , ^i(i))tG[o.i] is Markovian. The covariance function of 
the Gaussian process (X„(t))(£[o,i] is given by 



(s-u)"-i (t-u) 



n-l 



■ du. 



(n-l)! (n-l)! 

In order to apply Theorem 12.11 we decompose into the form ()2.ip . We have for, e.g., s < t, 



Cx„(s>*) = 



(s — uY' 



DP 



E 



L fc=o 



n — 1 



i-uy 



i—l — k^k 



du 



k=0 

n-l 



j-k /-s 



(s-u) 



n— 1 „,n — l — k 



{n~iy. (n-l-fc)! 

2n— 1 — /c 



du 



k=0 



(2n - 1 - A:)! kl 
We then obtain the following representation: 

with, for any k G {0, 1, . . . , n — 1}, 

= (-1)"-!-'= 



{2n-l-k)r = 



We state below a result of [10] that we revisit here by using Theorem 12. II 

Theorem 3.1 Let u be a fixed continuous function on [0, 1]. The function v defined on [0, 1] by 



v{t) = / c^^{s,t)u{s) ds 



is the solution of the boundary value problem 



'v(2») = (_i)«M on [0,1], 
i;(')(0) = for ie {0,1,..., n-l}, 

x;W(l) = /ori G {n,n + l,...,2n-l}. 



(3.1) 



Proof Let us check that the conditions ()2.2p and ()2.3p of Theorem 12. H are fulfilled. First, we have 

n-l 

j2[p,{t)^l:\t)-^\:\t)Mt)^ 



k=0 
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(-1)^ 



i-iy 



j2,n—l — i — k -jJi 



^{*^->'} (2n - 1 - fc)! (k - i)\ ^{'=<2"-i-»} (2n - 1 - I - fc)! A;! 



4-2n — 1 — 2 



(2n - 1 - i)! 

2ri-l-i 



k—i 

/ n— 1— i 



2?i - 1 - z 



k — i 



(2n-l-i)A(n-l) 



E (-1) 



fc=0 



2n - 1 - z 



A;=0 



, /2?i - 1 - z 



E(-i)* 



2n - 1 - z 



Performing the transformation A; i— > 2n — 1 — i — A; in the first sum lying within the last equality, 
we get 



n — l — i 



i+k 



em: 

fc=0 

2n-l-i 

= E (-1)* 

k—n 

and then 



. /2n - 1 - i 



-EM) 



fe=0 



2n - 1 - i 



2?i - 1 - i 
2?T. — 1 — « — A 



E(-i) 



. /2n - 1 - i 



fe=0 



2)1-1-1 



= E (-1)* 



2n - 1 - i 



— <5j,2n-l 



n—L 



fc=0 



On the other hand, setting 



v = {-iy 



we clearly see that, for any fee {0, 1, . . . , n — 1}, 

V^k = V^k = and l?^(^fc(0) =PiV^fc(l) =0 for i G {0, 1, . . . , n - 1}. 

Consequently, by Theorem 12.11 we see that the function v solves the boundary value problem p.ip . 
The uniqueness part will follow from a more general argument stated in the proof of Theorem 14.31 
□ 



4 Various bridges of the (n — l)-fold integral of Brownian motion 

In this section, we construct various bridges related to {Xn{t))tfz[o,i]- More precisely, we take 
{Xn{t))te[o,i] conditioned on the event that certain derivatives vanish at time 1. Let us recall that 
all the derivatives at time naturally vanish: X„(0) = X„_i(0) = • • • = ^^^2(0) — Xi{0) — 0. 

For any m G {0, 1, . . . , n}, let J = {ji, J2, ■ • • , im} be a subset of {1, 2, . . . , n} with 1 < ji < j2 < 
■ ■ ■ < jm < n and the convention that for m = 0, J = 0. We see that for each m € {0,1,..., n}, 
we can define (") subsets of indices J, and the total number of sets J is then V" n f") —2". Set 

\ni/ ' A^?n— \ni/ 

for any t G [0, 1] 

Y{t) = {x^{t)\x,ii) = 0, J e J) = (x„(i)|x,,(i) = x,,(i) = • • • = x,^{i) = 0). 

In this way, we define 2" processes (y(i))tg[o_i] that we shall call "bridges" of {Xn{t))te[o,i]- In 
particular, 

- for J = 0, we simply have (i^(t))te[oa] = (Xn{t))telo,i]; 

— for J = {1}, the corresponding process is the (n — l)-fold integral of Brownian bridge 

(X„(<)|Xi(l) = O),e[o,il = ( /* ^7^^d/3(.)) ; 

\Jq (n-i)'. /te[o,i] 
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— for J = {n}, the corresponding process is the "single" bridge of (n— l)-fold integral of Brownian 
Brownian: 



te[o,i] 



— for J = {1, 2, . . . , n}, the corresponding process is 

(X„(t)|X„(l) - X„_i(l) - • • . = = 0)4e[o,i]- 

This is the natural bridge related to the n-dimensional Markov process {X„{t), Xn-i{t), . . . , 
-'^i(*))te[o,i]- 

In this section, we exhibit several interesting properties of the various processes (^(i))te[o,i]- 
One of the main goals is to relate these bridges to additional boundary value conditions at 1. For 
this, we introduce the following subset / of {0, 1, . . . , 2n — 1}: 

I ={n- J)U [{n, n + 1, . . . , 2n - 1}\( J + n - I)] 



with 

n- J = {n- j,j e J} = {n~ ji,. .. ,n- j„J, 
J + n-l = {j + n-l,j e J} = {ji+n-1,... J„i +n-l}. 



The cardinality of / is n. Actually, the set / will be used further for enumerating the boundary 
value problems which can be related to the bridges labeled by J. Conversely, / yields J through 



J = (n-/)n{l,2,...,n}. 



In the table below, we give some examples of sets / and J. 



I 


{0,l,...,n-l} 


{n, n + 1, . . . , 2n — 1} 


{n - l,n + 1, . . . ,2n - 1} 


{0,n + l,...,2r7.-2} 


J 


{1,2,. 





{1} 


{n} 



4.1 Polynomial drift description 

Below, we provide a representation of (^(i))te[o,i] by means of (X„(i))jg[o ij subject to a random 
polynomial drift. 

Theorem 4.1 We have the distributional identity 

\ jgj /t6[0,l] 

where the functions Pj, j G J, are Hermite interpolation polynomials on [0, 1] characterized by 

Pf(0) = /ori G {0,l,...,n-l}, (4.1) 

^j'^(l) = '5j-„_, for tel. 

Remark 4-1 In the case where n = 2, we retrieve a result of [Bj. Moreover, the conditions (|4.ip 
characterize the polynomials Pj, j G J. We prove this fact in Lemma lA. II in the appendix. 
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Proof By invoking classical arguments of Gaussian processes theory, we have the distributional 
identity 

V /te[04] 

where the functions Pj, j G J, are such that E[y(t)Xfc(l)] = for all k £ J. We get the linear 
system 

^E[X,(l)Xfc(l)] P,it) = E[X„(t)Xfc(l)], keJ. (4.2) 

We plainly have 

(.?-!)! 



Then, the system (|4.2p writes 



f^^j + k-iij-iy. J, [n-iy. 

The matrix of this system (l/(j + fc — 1))^ is regular as it can be seen by introducing the related 
quadratic form which is definite positive. Indeed, for any real numbers Xj, j e J, we have 

. ^^^'^ =0 if and only if £ J, a;^ = 0. 



3,keJ " \j,feGJ 

and 



?■ + A: - 1 



Thus, the system (|4.2p has a unique solution. As a result, the are linear combinations of the 
functions t i-> — — u)'^^^ du which are polynomials of degree less than n + k. Hence, 

Pj is a polynomial of degree at most 2?! — 1. 

We now compute the derivatives of Pj at en 1. We have Pj*'(0) = for i G {0, 1, . . . , n — 1} 
since the functions t f^{t — u)^~^{l — u)^~^ du plainly enjoy this property. For checking that 
Pj^\l) = Sj^n-i for i £ I, we successively compute 
— for i e {0, 1, . . . , n — 1}, 

^E[X„(t)Xfc(l)] = E[Xn-^{t)Xk{l)]■, 



— for i = n — 1, 



— for i = n, 



^E[X„(i)X,(l)] = E[B{t)X,{l)] = ^ ^±^^L_du 



d* (1 - t)''^'^ 

— E[X„(i)X,(l)]- ^ > ■ 



dP ' ' " (fc- 1)! ' 
for i G {n, n + l,...,ri + fc — 1}, 

— E[X„(t)X,(l)] = (-1)'+" ^ . 

di* (k + n — z — Ij! 

^E[X„(t)Xfe(l)] = (-1)^-1; 



for I = fc + 7i — 1 , 



— for i > fc + n. 



^E[X„(t)Xfe(l)]=0. 
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Consequently, at time t — 1, for i>n, 

^E[X„(t)Xfc(l)][^^ = {-lf-H,,k+n-i. (4.3) 
Now, by differentiating (|4.2p . we get for i € {0, 1, . . . , n — 1}, 

Y,¥.[Xj{l)Xu{l)]Pf{l) = E[X„_,(l)Xfc(l)], fc e J. 

In particular, if i 6 (n — J), this can be rewritten as 

^E[X,(l)Xfc(l)] Pf\l) = ^,5,,„_,E[X,(l)Xfc(l)], fc e J, 
jeJ je,J 

which by identification yields P^^*''(l) = Sj^n-i- Similarly, for i G {n, n + 1, . . . , 2n — 1}, in view 
of (|4.3p . we have 

^E[X,(l)Xfc(l)] /^(^'(l) = (-l)^-i(5,,fc+„_i, keJ. 
In particular, if i G {n, n + 1, . . . , 2n — 1}\( J + n — 1), we have Si^k+n-i = for k £ J, and then 

^E[X,(1)X,(1)]P,«(1) = 0, fee J, 

which by identification yields Pj^\l) = = Sj^n-i- The proof of Theorem 14. II is finished. □ 
4.2 Covariance function 

Let Cy be the covariance function of (i^(i))te[o.i] ■ Cy.{s,t) = K[Y{s)Y{t)]. In the next theorem, we 
supply a representation of of the form (12.11) . 

Theorem 4.2 The covariance function of (i^(i))t6[o,i] admits the following representation: for 
s,iG[0,l], 

n-l 
fc=0 

with, for any k £ {0, 1, . . . , n — 1}, 

Moreover, the functions ijjk, k £ {0, 1, . . . , n — 1}, are Hermite interpolation polynomials such that 

4'k\0)=kk /or^e{0,l,...,n-l}, 
^V;i'^(l)=0 for I el. 
Proof We decompose y(t) into the difference 

Y{t) = Xnit) - Z{t) with ^ 

We have 

c^(s,i) - E[X„(s)X„(t)] +E[Z(s)Z(<)] - E[X„(s)Z(t)] - E[Z(s)X„(t)] 
= c,Js,0+ E IE[^,(l)^fe(l)]i',(s)Pfe(t) 

-Y,nXn{t)X,{l)] P,{S) -Y,nXn{s)X,{l)] P,{t). 
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By definition (I4.2p of tlie -Pj's, we observe tliat 

E[X,{l)Xk{l)] PAs)Pk{t) -Y,nXn{s)Xu{l)] Pk{t) 



keJ 



k£j VjSJ / 

Tlicn, we can simplify c-^(s,t) into 

(s, t) = (s, i) - E (O^j (1)] (5)- 

Since the covariance functions and c^^ are symmetric, we also have 

(s, = (s, - E E[X„(s)X, (1)] P,{t). 
je-r 

Let us introduce the symmetric polynomial 

Q{s,t) = EE[X„(s)X,(1)] Pj{t) =Y,nXn{t)X,il)] P,{s). 

It can be expressed by means of the functions ipk, V'fe's as follows: 

ji-i 

Q{s, = E E A t)4-'^ (1)P, ( W t). 



We can rewrite Cy (s, t) as 
and then, for s,t G [0, 1], 



fc=0 



Cv(s,0 = E'^fc(sAt)Vife(sVt) with i)k{t)^Mt)-J2^k i^)P3it)- 



k=0 



We immediately see that 'ipk is a polynomial of degree less than 2n such that V'[*''(0) = Si^k for 
is {0, 1, . . . , n — 1} and, since Pj*^(l) = 5j^n-i, 



We deduce that 



if i £ in — J), 



= Vfe \l) =0 if i e {n, n + 1, . . . , 2n - 1}\(J + n - 1). 
Then '^[^''(l) = for any i £ I. This ends up the proof of Theorem 14.21 □ 
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4.3 Boundary value problem 

In this section, we write out the natural boundary value problem which is associated with the 
process (i^(i))te[o,i] • The following statement is the main connection between the different boundary 
value conditions associated with the operator d?" /dt^" and the different bridges of the process 
(X„)jg[o,i] introduced in this work. 

Theorem 4.3 Let u he a fixed continuous Junction on [0, 1]. The function v defined on [0, 1] by 

v{t) = / (s, t)u(s) ds 

is the solution of the boundary value problem 

w(2«) ^ on [0,1], 

uW(0)=0 /or i e {0, l,...,n- 1}, (4.4) 

uW(l)=0 fori el. 

Proof • First step. Recall that (s, t) — (s, t) — Q{s, t). We decompose the function v into the 
difference w — z where, for t 6 [0, 1], 

w{t) = / Cj^^{s,t)u{s) ds, z{t) — / Q{s,t)u{s) ds. 
Jo " Jo 

We know from Theorem O that w^^") ^ wW(0) = ^{s,0) = for i G {0, 1, . . . , n - 1} 

and ly*^*' (1) = for i G {n, n + 1, . . . , 2n — 1}. Moreover, the function 1 1— Q{s, t) being a polynomial 
of degree less than 2n, the function z is also a polynomial of degree less than 2n. Then z^^n) = q, 
(0) = for i e {0, 1, . . . , n - 1} and 

^(2«) ^ ^(2«) _ ^(2n) ^ „W(o) = yj(i){Q) _ ^W(o) =0 for i £ {0, 1, . . . , n - 1}. 

On the other hand, we learn from (|2.5p that, for i G {0, 1, . . . , 2n — 1}, 

and then, for i G {0, 1, . . . , n — 1}, 

z(;«(l)= / E[X„Gs)X„_,(l)]u(s)ds. 
Jo 

We also have, for any i e {0, 1, . . . , 2n — 1}, 

^(s,l) = ^E[X„(s)X,(l)]P,«(l) = l{,g(„_^)jE[X„(s)X„_,(l)]. 

As a result, we see that 

^^'Hl) = ^(s, = (1). 

This implies that for i € {n - J), z(*)(l) w^\l) and for i e {n, n + 1, . . . , 2n - l}\(J + ri - 1), 
z(*)(l) = = Then = for any i £ /, that is = 0. The function w is a 

solution of (gH]). 



12 



• Second step. We now check the uniqueness of the solution of (|4.4p . Let vi and V2 be two solutions 
of Vv = u. Then the function w = vi — V2 satisfies Vw — 0, w*^*'(0) = for i G {0, 1, . . . , n — 1} 
and u)(*^(l) = for « e /. We compute the following "energy" integral: 



1 

u;(")(i)2dt = (-1) 



-1 



+ (-1)" / w(t)u;('")(t)di 







1=0 



We have constructed the set / in order to have w*^*^ (l)w(^" ^ (1) = for any i e {0, 1, . . . , n — 1}: 
when we pick an index i G {0, 1, . . . , n — 1}, either i G / or 2n — 1 — i G /. Indeed, 

- if i G / n{0,l,...,n- 1}, = 0; 

- if i G {0, 1, . . . , n — by observing that {0, 1, . . . , n — = {0, 1, . . . , n — ^ J), we 
have 2n - 1 - i G + 1, . . . , 2n - 1}\( J+ n - 1). Since {n,n + l, . . . , 2n - 1}\( J + n - 1) C /, 
we see that 2n — 1 — i G / and then it;(2"-i-0 — 

Next, w(")(t)2dt = which entails w;^") = 0, that is, It; is a polynomial of degree less than n. 
Moreover, with the boundary value conditions at 0, we obtain w = or wi = V2. 
The proof of Theorem 14.31 is finished. □ 



Remark ^.2 We have seen in the above proof that uniqueness is assured as soon as the boundary 
value conditions at 1 satisfy w'*-' (1) = for any i G {0, 1, . . . , n — 1}. These conditions 
are fulfilled when the set / = {ii, 12, . . . , i„} C {0, 1, . . . , 2n — 1} is such that ii G {0, 2n — 1}, 
12 € {l,2n — 2}, 13 G {2,2n — 3}, i„ G {n — This is equivalent to say that / and 

(2n — 1 — /) make up a partition of {0, 1, . . . , 2n — 1}, or 

2n - 1 - / = {0, l,...,2n- 

In this manner, we get 2" different boundary value problems which correspond to the 2" different 
bridges we have constructed. We shall see in Section [6] that the above identity concerning the 
differentiating set / characterizes the possibility for the Green function of the boundary value 
problem (|4.4p to be the covariance of a Gaussian process. 



4.4 Prediction 

Now, we tackle the problem of the prediction for the process (i^(i))te[o.i] ■ 

Theorem 4.4 Fix to G [0, 1]. The shifted process {Y(t + io))te[o,i-to] (admits the following repre- 
sentation: 

{Y{t + to))te[o.i-to] - [Yto{t) + J2Q^>to{t)Y^'Hto)] 

V «=0 / tG[0,l-to] 

where 

Ytoit) = Xjt) ~J2P^,,Jt)X,{l - to). 

The process {Xn{t)) is a copy of {Xn{t))te[o,i-to] which is independent of {Xn{t))te[o,to]- 
The functions Pj.to, j G J , Qi.to, i G {0, 1, . . . , n—l}, are Hermite interpolation polynomials 



on [0,1 — to] characterized by 



^7".tn — 



i,*o 

i^S(0) = /ortG{0,l,...,n-l}, 
j(0 



VPll,{\-to) = k,n-, for Lei, 
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fO^^") = 
Qmo(0) = ^m /orie{0,l,...,n-l}, 

Actually, these functions can be expressed by means of the functions Pj, j G J, and ipi, i G 
{0, 1, . . . ,n — 1}, as follows: 



In other words, the process (^to(^))igjQ 2_t ] bridge of length (1 — to) which is independent of 



Proof • First step. Fix to G [0, 1]. We have the weU-known decomposition, based on the classical pre- 
diction property of Brownian motion stipulating that (i?(t-|-to))te[o,i-to] — (^to(0+^(^o))tg[o i-to] 
where (^tg (t))^gjj^ is a Brownian motion independent of (B(t))tg[o_f(,] , 



(X„(t + to))fe[o,l-to] = ( ^n{t) + X/ ^ Xn-k{ti) 



k=0 



te[o,i-to] 



with Xn{t) — Jq ^^(^n-liy. d-B(s). Differentiating this equality {n — j) times, j S {1,2,..., n}, we 
obtain 

X,{t + t^)^Xj{t)+ J2 7--l—X\ ^n-k{to) 



k—n—j 



{j + k-n)\ 



Therefore, 

r(t + to) 

= X„(t + to) -^F,(t + to)X,(l) 



k=0 



{j + k — n)! 



Xa{t) -Y,Pjit + to)Xj{l-to) 



k=0 L jeJ:j>n-k ' 



Xn-k{to)- 

(4.5) 



We are going to express the X„_fc(to), fee {1, . . . , n}, by means of the Xj{l — to), 6 J. We have, 
by differentiating ()4.5p i times, for i G {0, 1, . . . , n — 1}, 



r«(^ + io) 



X„_,(t) - ^ Pf\t + to)l,(l - to) 

fe=0 L ^ ' jGJ:j>n-fc ' 



Xn-k{to)- 



For t = 0, this yields 



Y(^\to) = -J2pj'\to)X,{l-t,) + J2 
jeJ k=o 



jGJ:j>n-k 



{j + k — n)l 



Xn-k{to). (4.6) 
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Set for i,fc G {0,l,...,n- 1} 

and let us introduce the matrix 
together with its inverse matrix 



{j + k-ny. ' 



A — {aik)o<i,k<n- 



B — A ^ — {bik)o<i,k<n-l- 

The equahties ()4.6p for i e {0, 1, . . . , n — 1} read as a hnear system of {n — 1) equations and (n — 1) 
unknowns: 



a,kXn-kito) = Y'^'\t„) + Pjito)X,il - to), I e {0, 1, . . . , n - 1}, 



which can be rewritten into a matrix form as 



A 



Xn-l{to) 
X2{to) 



The solution is given by 

/ Xnih) \ 
Xn-l{tQ) 



^2 (to) 



B 



r("-2)(to) 
Vr(«-i)(to)y 



Y^jil-to)B 



p;{to) 



\Pt''\to)J 



( Pjito) \ 

p'Ah) 

Pj"-'\to) 
\Pt~'\to)J 



and we see that Xn-k{to), k e {0, 1, . . . , n — 1}, is of the form 



1=0 



i=0 



1,(1 -to). 



Therefore, by plugging ()4.7p into ()4.5p . we obtain 

Y{t + to) = Xn{t) - Y Pjii + ^o)^,(l - to) 

J 6./ 

n-1 



/c=0 
n-1 



jeJ:j>n-k ' 



^ 5..,y«(to) + E E &fe«^^''(^o) ^,(1 - to) 

i=0 jeJ \i=0 / 



Fj(t + io)- E ^^^Pf^to) 



a<i,k<n-\ 



= x,,it)-j2 

(t^ (l-to)^+^-" ^ 



fc! ^ (7 + fc-n)! 



1,(1 -to) 



+ E 

0<j,fc<n-l 



ik 

E 



(1 - to)^+'^"" 



k\ ^ (j + k-n)\ 



Pj{t + to) 



(4.7) 



yW(to). 
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Finally, Y{t + to) can be written as 

n-l 

Y{t + to) = Yt, it) + J2 Q^,to it)Y^'^ (to) 



i=0 



where 



and 



with 



fe=0 



y 



(1 - ioP+^-" 



fc! . /-t , (j + fc-n)! 



Pj{t + to) 



Yt,{t) = 1„(0 - ^ P,,t„ - to) 



0<!,A:<n-l 



je.J:]>n-k ' 



^P,it + to)-J2pj'\to)Q^,toit) 



i=0 



• Second step. We easily see that the functions Pj^to and Qi^to are polynomials of degree less than 
2n. Let us compute now their derivatives at and ^o- First, concerning Qij^ we have 



j-k — L 



(k-i) 



■1 



{k>i.} 



jeJ:]>n-k 



{3 + k — n)\ 



Choosing t — and recalling the definition of a^k and the fact that the matrices {aik)o<i,k<n and 
(&ife)o<i,fc<n are inverse, this gives for t G {0, 1, . . . , n - 1}, 



j^J'j^n—k 



{j + k- n) 



Choosing i = 1 — io, we have for t G / 



i;=0 



4fe>t} 



j&J-j>n — k 



U + k-ny. 



(4.8) 



By Theorem 14.11 '^e know that P'-'^-'(l) = (5j.„-t for t G /; then 



E 



jeJ:j>n-k 



, (j + fc-?i)! 



_ ,V l{fc>4l{te(n-J)}- 



Observing that, ifi<fc<n — 1, the conditions i G (n — J) and t G / are equivalent, we simply 
have 

(1 - to)^-+^^" _ (1 - to)^- 



E 



jeJ:j>n-k 

which immediately entails, by 



Next, concerning Pj''tg, we have 



(i + - n)\ 



P^>(1) 



(fc-.)! 



1 



{fc>t} 



(1 - io) = for L G /. 



i=0 
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Choosing t = 0, this gives for t e {0, 1, . . . , n — 1}, since Q-'j^(O) = S^^i, 

n-1 

P^W (0) = p^^\to) - J2 Pj'\to)Q[i{0) = 0. 

1=0 

Choosing t = 1 — to, we have for t G /, since pj'^'(l) = 5,„„_j and Qi'f(,(l — to) = 0, 

(1 - to) = Pj'\l) - E ^^'H*0)Q!:i(l - to) = 5.,n-,. 



The polynomials Pj^to (resp. Qi,to) enjoy the same properties than the Pj's (resp. the V'i's)! re- 
garding the successive derivatives, they can be deduced from these latter by a rescaling according 
as 

P,,i„(t) = (1 - io)«-^P,. (^^^ (resp. g,,i„ = (l-io)V.(^^)). 

It is then easy to extract the identity in distribution below, by using the property of Gaussian 
conditioning: 

[ X„(t) -^P,,t„(t)l,(l - to)) - (In (t) 1 1,(1 - to) = 0,J e J)t^[o,l-toy 



te[04-to] 

Theorem 14. 41 is established. □ 



5 Example: bridges of integrated Brownian motion (n — 2) 

Here, we have a look on the particular case where n = 2 for which the corresponding process 
(X„(t))jg[o,i] is nothing but integrated Brownian motion (the so-called Langevin process): 



X2{t) = / B{s)ds. 
Jo 



The underlying Markov process is the so-called Kolmogorov diffusion {X2{t) , Xi{t))t^[f).i]- AH the 
associated conditioned processes that will be constructed are related to the equation i;^"') (t) — u{t) 
with boundary value conditions at time 0: v{0) = v'{0) = 0. There are four such processes: 

— (X2(t))jg[o,i] (integrated Brownian motion) ; 

— (X2(t)|Xi(l) = 0)tg[o,i] (integrated Brownian bridge); 

— {X2{t)\X2{l) — 0)tg[o,i] (bridge of integrated Brownian motion); 

— {X2{t)\Xi{l) — ^2(1) = O)tg[o,i] (another bridge of integrated Brownian motion). 

On the other hand, when adding two boundary value conditions at time 1 to the foregoing equation, 
we find six boundary value problems: v{l) = v'{l) ~ 0, t;(l) — v"{l) — 0, v{l) — v'"{l) — 0, 
v'{l) = v"{l) = 0, w'(l) = v"'il) = 0, v"{l) = v"'{l) = 0. Actually, only four of them can be 
related to some Gaussian processes-the above listed processes-in the sense of our work whereas 
two others cannot be. 

For each process, we provide the covariance function, the representation by means of integrated 
Brownian motion subject to a random polynomial drift, the related boundary value conditions at 
1 and the decomposition related to the prediction problem. Since the computations are straight- 
forward, we shall omit them and we only report here the results. 

For an account on integrated Brownian motion in relation with the present work, we refer the 
reader to, e.g., [5] and references therein. 
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5.1 Integrated Brownian motion 

The process corresponding to the set J = is nothing but integrated Brownian motion: 



(^2(i))te[o,il = (£^(s)ds 



te[o,ij 



The covariance function is explicitly given by 



c{s,t) = ^[sAtf [3(s Vt) -sM]. 

This process is related to the boundary value conditions at 1 (/ = {2,3}): v"{l) = v"'{l) = 0. The 
prediction property can be stated as follows: 

{X2{t + to))te[o,i-to] = {Mt) + X2{to) +tX,{to))^^^^,_^^y 
5.2 Integrated Brownian bridge 

The process corresponding to the set J = {1} is integrated Brownian bridge: 
{y{t))te[o,i] = ( fB{s)ds B{1) = o) =( fmds) = {X2{t)\X,{l) = O)te[o,iv 

\J0 /f6[0,l] V-^O /t6[0,ll 

This process can be represented as 



iY{t)),^^o,i] = (x2{t)-lt'X,{l)) 

V ^ / te[o,i] 



The covariance function is explicitly given by 

c(s, t) = - [s A [3(s Vt)-sAt]-- s^t^. 
6 4 

The process (X{t))te[o,i\ is related to the boundary value conditions at 1 (/ = {1,3}): v'{l) 
v"'{l) = 0. The prediction property says that 

+2 



te[o,i-to] 



5.3 Bridge of integrated Brownian motion 

The process corresponding to the set J = {2} is the bridge of integrated Brownian motion: 

{Y{t))t^lo,i] = f / Bis) ds f Bis) ds^o] = (X2(f)|X2(l) = O)te[o,i]. 

V^o Jo /te[o,i] 

The bridge is understood as the process is pinned at its extremities: 5^(0) = 5^(1) = 0. This process 
can be represented as 

iYit)),^,o,i] = (x2it)-lt'i3-t)X2il) 



te[o,i] = I ^2yt) - 

V /te[o,i] 

The covariance function is explicitly given by 

c(s,t) = i[sAi]2[3(.sVt) - sAt]- ^ .5^*2(3 - s)(3 - t). 
6 12 

The process (5^(i))te[o,i] is related to the boundary value conditions at 1 (/ = {0,2}): vil) - 
v"il) =0. The prediction property says that 

(y(i + to))t6[o,i-to] 

(y m , - 3(1 -to)^^ + 2(1-^0)3 t^-3il-to)t-' + 2il-to)H \ 
= V''^'^ ^ "^^''^ ^ ^'°V.e[o,i-M 
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5.4 Other bridge of integrated Brownian motion 



The process corresponding to the set J = {1, 2} is another bridge of integrated Brownian motion 
(actually of the two-dimensional Kolmogorov diffusion) : 



iY{t))te[o,i] = { 1^ B{s)ds 



B{s)ds = B{1) =0 = {X2{t)\X2{l) = = O),e[o,i]- 

/ te[o,i] 



The bridge here is understood as the process is pinned at its extremities together with its deriva- 
tives: Y{0) = Y'{0) = Y{1) = Y'{1) = 0. This process can be represented as 



{Y{t))te[o,i] = {X2{t) - t\t - l)Xi(l) - t\3 - 2t)X2(l)),^[o_,j . 
The covariance function is explicitly given by 

c{s,t) = ^[sAt]'^ [1 - s Vi]2 [3(s Vf) -sAt-2st]. 



The process (i^(i))te[o,i] is related to the boundary value conditions at 1 (7 = {0,1}): v{l) = 
t;'(l) = 0. The prediction property says that 

[r[t + to))t£[o,i_t„] - I Yto [t) + T- .3 y{to)-\ 7Z y [to) ] 



5.5 Two counterexamples 

• The solution of the problem associated with the boundary value conditions v{l) = v"'{l) = 
(which corresponds to the set Ii = {0, 3}) is given by 

v{t)= [ Gi{s,t)u{s)ds, ie[0,l], 

where 

Gi{s,t) = ^[sAt]'^ [3(s Vt) -sAt] + ^s^t^(s- 3). 
6 6 

• The solution of the problem associated with the boundary value conditions v'{l) = v"{l) = 
(which corresponds to the set I2 = {1, 2}), is given by 

v{t)= [ G2{s,t)u{s)ds, ie[0,l], 
Jo 

where 

G2(s,t) = llsAt]^ \3(sVt) - s At] + lsh^(t - 3). 
6 6 

Wc can observe the relationships Gi{s,t) = G2(t, .s) and I2 = {0, 1, 2, 3}\(3 — Ii). The Green 
functions Gi and G2 are not symmetric, so they cannot be viewed as the covariance functions of 
any Gaussian process. In the next section, we give an explanation of these observations. 
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6 General boundary value conditions 

In this last part, we address the problem of relating the general boundary value problem 

1,(2") ^ on [0,1], 

v{0) = w'(0) = • • ■ = (0) = 0, (6.1) 

= = , . . = vM(l) ^ 0, 

for any indices ii, J2, . . . , in such that < ii < 12 < • • • < < 2n — 1, to some possible Gaussian 
process. Set / = {ii, 12, . . . , i„}. We have noticed in Theorem 14.31 and Remark 14.21 that, when 
I satisfies the relationship 2n — 1 — / = {0, 1, . . . , 2n — 1}\^, the system ()6.ip admits a unique 
solution. We proved this fact by computing an energy integral. Actually, this fact holds for any set 
of indices /; see Lemma [A. II 

Our aim is to characterize the set of indices / for which the Green function of (j6.ip can be 
viewed as the covariance function of a Gaussian process. A necessary condition for a function of 
two variables to be the covariance function of a Gaussian process is that it must be symmetric. So, 
we shall characterize the set of indices / for which the Green function of ()6.ip is symmetric and we 
shall see that in this case this function is a covariance function. 



6.1 Representation of the solution 



We first write out a representation for the Green function of (j6.ip . 

Theorem 6.1 The boundary value problem i6.1\) has a unique solution. The corresponding Green 
function admits the following representation, for s,t E [0, 1]; 



GAs,t) = (-iri^,<,jiL_f)__ _ 



(2n- 



(2n-l-i)! 

LEI ^ ' 



where the Ri^u, l E I , are Hermite interpolation polynomials satisfying 



(6.2) 



i?S(0) = /orie{0,l,...,n-l}, 
i?^:l(l) = '^M forieL 

characterize the polynomials l £ L We prove this fact in 



Remark 6.1 The conditions 
Lemma lA. II in the appendix. 

Proof Let us introduce the functions vi and V2 defined, for any t G [0, 1], by 



v,{t)^{-iY 



[t-sf 



(2n-l)! 



i(s)ds and V2{t) — v{t) — vi(t). 



We plainly have 7;^"^ = (— l)"u and wi(0) = v[{(S) = ■ ■ 
V solves the system (16.11) if and only if the function V2 satisfies 



(0) = 0. Therefore, the function 



„(2n) _ 



on [0,1], 



««(0)=0 



4')(1) = (-1)"-! [' ^L^)^lllu{s) 



for i G {0,l,...,n- 1}, 
ds for i Cz I. 



(6.3) 



Referring to Lemma lA. II the conditions (16.31) mean that V2 is a Hermite interpolation polynomial 
which can be expressed as a linear combination of the i G /, defined in Theorem l6.1l as follows: 



V2 



{t)^j2^i'\i)R,At) = {-ir~' f 



(l-,s) 



^ (2n 



2ri-l-i. 



■Ri.At) 



u{s) ds. 
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Consequently, the boundary value problem (|6.1[) admits a unique solution v whicli writes 

»i 



v{t) = vi{t) + V2{t) = G,{s,t)u{s)ds, te[0,l], 
Jo 

where Gj{s,t) is defined in Theorem 16. II The proof is finished. □ 

We now state two intermediate results which will be used in the proof of Theorem [ 

Proposition 6.1 Let Ii and I2 be two subsets of {0, 1, . . . , 2n — 1} with cardinality n. If the sets 
Il and I2 are different, then the corresponding Green functions Gj^ and Gj^ are different. 

Proof Suppose that Ii 7^/2- If we had G^^ — G^^ , for any continuous function u, the function v 
defined on [0, 1] by 

v{t)^ / G, {s,t)u{s)ds^ / G, {s,t)u{s)ds 



would be a solution of the equation u^^") — (— satisfying v^^\0) =0 for z G {0, 1, . . . , n — 1} 
and = ioi i e h U h- Since h ^ I2 and since h , I2 have the same cardinality, there 

exists an index io which belongs to l2\Ii- For this io, we would have {d'^"Gj_^/dt^"){s, 1) = for all 
s G (0, 1), or equivalently, 

(1 - ,)2n-l-t ^ (1 - 5)^"'!-'" 

4- {2n-l-L)\ ""^-'^'^ (2n-l~zo)! ' 
ttii 

This is impossible since the exponent (2n — 1 — iq) does not appear in the polynomial on the 
left-hand side of the foregoing equality. As a result, Gj^ ^ G,^ . □ 

Proposition 6.2 Let Ii and I2 be two subsets of {0, 1, . . . , 2n — 1} with cardinality n. The re- 
lationship Gj_^{s,t) — G^^{t,s) holds for any s,t £ [0,1] (in other words, the integral opera- 
tors with kernels G^^ and Gj^ are dual) if and only if the sets Ii and I2 are linked by I2 = 
{0,l,...,2n-l}\(2A-l-/i'). 

Proof We have, for any s, t G [0, 1], 

(-1)"[G,(M) = - E ^n'-l-.)! ^---^^^ + ^ 2n- l-.)! ^--^^^- 

Set, for any s, i £ [0, 1], 

S{s,t) = {-ir[G,^{s,t) - G,^{t,s)]. 
The polynomial S has a degree less than 2n with respect to each variable s and t and satisfy 

Q^l+^2S (i_s)2n-l-»l-»2 
(s,i) = (-ir%+.2<2„-l} ^ ' 



..^ , fl _ „^2rl-l-t-^2 

(i-i)^ 

^{ii<2n-l-t} 

te/2 

In particular, 

— for t = 1 and ii G Ii, 12 G {0, 1, . . . , 2n — 1}, 



2n — 1 — t — ti 



(s,l)-(-l)^^%+,,<2„-l} ^ ^ 



E(l — s)2n-l-t-t:2 
l{t2<2ri-l-t} 775 ^ -T7 (51,41 
^ (2n — 1 — i — 22)' 

te/i 

te/2 

= (-l)^^%e(2„-i-/2)}i?£tGs); (6.4) 
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— for {s,t) = (1,0) and ii G {0, 1, . . . , n - 1}, 12 £ h, 

())-( ^Y^+^ l{ti+»2<2»-l} , iV'iV %<2n-l-4 „ _ , . 

as'^at.i ^1'"^-^ 1^ (2n-l-zi-z2)! +^ ^> 2_.(2n-l-.-zi)! ' ^ ^ ^ 

— for (s, t) = (0, 0) and ii, 12 £ {0, 1, . . . , n — 1}, 

P)il+i2 Q 

^-^(0,0) = (-ir5,+.„2„-i=0. (6.6) 

Now we are able to establish the statement of Proposition 16.21 

— If /2 7^ {0, 1, . . . , 2n — l}\(2n — 1 — /i), there exists ii e Ii such that ii G (2n — 1 — 12) and 
then, in view of ()6.4p . 

The polynomial S* cannot be null, that is, there exist s, t G [0, 1] such that G^^ (s, i) 7^ G,^ {t, s). 

— If /2 = {0, 1, . . . , 2n — l}\(2n — 1 — Ii), for any ii G /i, we have ii ^ (2n — 1 — 12) and then, in 
view of ()6.4p . 

|^(s,l)=0 foriiG/i. 

Put, for any ii G {0, 1, . . . , n — 1}, Si^ (s) = -^—i - (s, 0). The polynomial S'i^ has a degree less 
than 2n. By (|6.5p and (|6.6p . we have 

^,^;^)(0)=0 forf2 e{0,l,...,n-l}, 
^S^l'\l)=0 for^2G/2, 

from which we deduce, invoking Lemma I A. 11 that all the polynomials Si-^ , ii G {0, 1, . . . , n — 1}, 
are null. Finally, the polynomial S has a degree less than to 2?! with respect to t and satisfies 

^^'^{s,0) = for ii G {0,l,...,n- 1}, 



(s, 1) = for ii G h. 



We can assert, by Lemma [A.li that S is the null-polynomial. 
The proof of Proposition 16.21 is finished. □ 

A necessary condition for Gj to be the covariance function of a Gaussian process is that it 
must be symmetric: Gj{s,t) — Gj{t,s) for any s,t G [0, 1]. The theorem below asserts that if the 
set of indices / is not of the form displayed in the preamble of Section HI that is / 7^ {0,1,..., 
2n — l}\(2n — 1 — /), the Green function Gj is not symmetric and consequently this function can 
not be viewed as a covariance function, that is we can not relate the boundary value problem (|6.ip 
to any Gaussian process. 

Theorem 6.2 The Green function G^ is symmetric (and it corresponds to a covariance function) 
if and only if the set of indices I satisfies 2n — 1 — / = {0, 1, . . . , 2n — 1}\/. 

Proof Set /' = {0, 1, . . . , 2n — l}\(2n — 1 — /). By Proposition 16.21 we see that Gj is symmetric if 
and only if G^, (s, t) = Gj (s, t) for any s, i G [0, 1], that is, by Proposition 16.11 if and only if / = /'. 
□ 



We made several verifications with the aid of Maple. Below is the program we wrote for this. 
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[> Green_f unction :=proc(n,setI) local V,M,S,T,P,setIcomp; 
V: =i->vector (n, [seq(binomial(k, i) ,k=n. . 2*n-l)] ) ; 
M: =llste->stackmatrix(seq(V(i) , i=liste) ) ; 

S:=(s,liste)->Matrix(l,n, [seq( (1-s) " (2*n-l-i)/(i ! * (2*n-l-i) !) ,i=liste)]) ; 
T:=t->Matrix(n,l, [seq([t"k] ,k=n. .2*n-l)]); 

P:=(s,t,liste)->multiply(S(s,liste) , inverse (transpose (M(liste) ) ) ,T(t)) [1,1] ; 
setlcomp: = [op({seq(i,i=0. .2*n-l)} minus {seq(2*n-l-i,i=setl)})] ; 
printCvalue of n'=n, 'differentiating indices set I_l'=setl); 
print ('Green function for s<t: GI_l(s,t)' 

=sort(simplify((t-s)*(2*n-l)/((2*n-l) !)-P(s,t,setI)) , [s,t] ,plex)) ; 
print ('Green function for s>t: GI_l(s,t)' 

=sort(simplify(-P(t,s,setI)) , [s,t] ,plex)) ; 
print (' symmetry test: GI_l(s,t)=GI_l(t,s)?' 

=evalb (simplify ( (t-s) * (2*n-l) /( (2*n-l) !) -P (s , t , setl) +P (t , s , set I) ) =0) ) ; 
print ('complementary set of 2n+l-I_l: I_2'=setIcomp) ; 

print ( 'difference between the two Green functions for s<t: GI_l(s,t)-GI_2(s,t) ' 

=sort (simplif y (-P(s , t , setl)+P(s , t , setlcomp) ) , [s , t] ,plex) ) ; 
print ( 'difference between the two Green functions for s>t : GI_l(s,t)-GI_2(s,t) ' 

=sort (simplif y(-P(s ,t , setlcomp) +P (s , t , setl) ) , [s , t] ,plex) ) ; 
print (' equality test between the two Green functions for s<t : GI_1 (s , t)=GI_2(t , s)? ' 

=evalb(slmpllfy((t-s)"(2*n-l)/((2*n-l) !) -P (s ,t , setl) +P(t , s , setlcomp) ) =0) ) ; 
print (' equality test between the two Green functions for s>t: GI_l(s ,t)=GI_2(t ,s)? ' 

=evalb (simplif y ( (t-s) * (2*n-l) / ( (2*n-l) ! ) -P (s , t , setlcomp) +P (t , s , setl) ) =0) ) ; 
end proc; 

To obtain the two Green functions associated with the sets 7i and I2 , Gj^ and G^^ , together 
with the cquahty test between them, run the command [> Green_function(n,i_i) ;. For instance, the 

return of the command [> Green_function(5, [2,3,5,6,8] ) ; is 

value of n = 3, differentiating indices set I_l = [1,4,5] 
Green function for s<t: GI_l(s,t) = 

-1/120 s-5 - 1/72 s-4 t-3 + 1/24 s-4 t + 1/18 s-3 t-3 - 1/12 3*3 t-2 
Green function for s>t: GI_l(s,t) = 

-1/120 s-5 + 1/24 s-4 t - 1/72 3*3 t-4 + 1/18 s-3 t-3 - 1/12 3-3 t-2 
symmetry test: GI_l(s,t)=GI_l(t,s)? = false 
complementary set of 2n+l-I_l: I_2 = [2,3,5] 
difference between the two Green functions for s<t: 

GI_l(s,t)-GI_2(s,t) = -1/72 s-4 t-4 + 1/72 s-3 t-4 
difference between the two Green functions for s>t: 

GI_l(s,t)-GI_2(s,t) = 1/72 s-4 fS - 1/72 s'S ^4 
equality test between the two Green functions for s<t : Gl_l(s,t)=Gl_2(t,s)? = true 
equality test between the two Green functions for s>t: Gl_l(s,t)=Gl_2(t,s)? = true 



6.2 Example: bridges of twice integrated Brownian motion (n = 3) 

Here, we have a look on the particular case where n = 3 for which the corresponding process 
{Xn{t))te[o,i] is the twice integrated Brownian motion: 

X^it) = s)B{s) ds = ^ (^jj B{si) dsi^ ds2. 

All the associated conditioned processes that can be constructed are related to the equation 
v'-^\t) = —u{t) with boundary value conditions at time 0: v{0) = v'{0) = v"{0) = 0. There 
are 2^ = 8 such processes. Since the computations are tedious and the explicit results are cumber- 
some, we only report the correspondance between bridges and boundary value conditions at time 1 
through the sets of indices / and J. These are written in the table below. 



conditioning set J 





{1} 


{2} 


{3} 


{1,2} 


{1,3} 


{2,3} 


{1,2,3} 


differentiating set / 


{3,4,5} 


{2,4,5} 


{1,3,5} 


{0,3,4} 


{1,2,5} 


{0,2,4} 


{0,1,3} 


{0,1,2} 



The Green functions related to the other sets cannot be related to some Gaussian processes. The 
sets are written in the table below with the correspondance I2 = {0, 1, 2, 3, 4, 5}\(5 — /i). 



differentiating set Ii 


{0,1,4} 


{0,1,5} 


{0,2,5} 


{0,3,5} 


{0,4,5} 


{1,4,5} 


differentiating set I2 


{0,2,3} 


{1,2,3} 


{1,2,4} 


{1,3,4} 


{2,3,4} 


{2,3,5} 
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6.3 Example: bridges of thrice integrated Brownian motion (n ~ 4) 

For n — A, only the 2'* = 16 following differentiating sets can be related to bridges: 

{0,1,2,3}, {0,1,2,4}, {0,1,3,5}, {0,1,4,5}, {0,2,3,6}, {0,2,4,6}, {0,3,5,6}, {0,4,5,6}, 
{1,2,3,7}, {1,2,4,7}, {1,3,5,7}, {1,4,5,7}, {2,3,6,7}, {2,4,6,7}, {3,5,6,7}, {4,5,6,7}. 



Appendix: Hermite interpolation polynomials 

.,n — 1}, and hi, i ^ I, be real numbers. There exists a unique 



Lemma A.l Let a^, i G {0, 1, , 

polynomial P such that 



p{2n) ^ 0^ 

pW(0)=a, /or {0,1, 
PW(1)=&, fori el. 



1}, 



(A.l) 



Remark A.l The conditions (lA.ll) characterize the Hermite interpolation polynomial at points 
and 1 with given values of the successive derivatives at up to order n — 1 and given values of the 
derivatives at 1 with selected orders in /. When / 7^ {0, 1, . . . , n — 1}, these polynomials differ from 
the usual Hermite interpolation polynomials which involve the successive derivatives at certain 
points progressively from zero order up to certain orders. 

Proof We look for polynomials P in the form P{t) — Y^Y=o^ ""J have 

j=o ^-^ ' 

We shall adopt the convention l/[(j — for i > j. The conditions ()A.ip yield the linear 

system (with the convention that i and j denote respectively the raw and column indices) 



Ci — a. 

I 2n-l 

V 



^ (j-z)! 

i=o ' 



if i £ {0,1, 
if i e /. 



1}, 



The (2n) x (2n) matrix of this system writes 



/ ((^ij) ie{o, 



A = 



3e{o,i, 



-1} 
-1} 



iei 

je{o.i,. 



I (^i,j) ie{o,i... 

3e{o,i,. 



(0).e{o.i...... 



iei 



Proving the statement of Lemma I A. 1 1 is equivalent to proving that the matrix A is regular. In view 
of the form of A as a bloc-matrix, we see, since the north-west bloc is the unit matrix and the 
north-east bloc is the null matrix, that this is equivalent to proving that the south-east bloc of A 



is regular. Let us call this latter Ao and label its columns C'j°\ j € {0, 1, . . . , n 

1 



1}: 



{n + j -i)! 



i€I 

je{o.i,. 



.(0) 



For proving that Aq is regular, we factorize Aq into the product of two regular triangular matrices. 
The method consists in performing several transformations on the columns of Aq which do not 
affect its rank. We provide in this way an algorithm leading to a LU-factorization of Aq where L 
is a lower triangular matrix and U is an upper triangular matrix with no vanishing diagonal term. 
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We begin by performing the transformation 



fcf 



c 



(0) 



ifi = 0, 



if je{l,2,...,n-l}. 



> n + j — ii 

The generic term of the column cj^\ for j G {1,2, . . . ,n — 1}, is 

1 1 1 i-ii 



1 



(n + j — z)! n + j — ii (n + j — z — 1)! n + j — ii (n + j — i)\ 
This transformation supplies a matrix Ai with columns C^-^-*, j G {0, 1, ... ,n — 1}, which writes 



{n-i)\ J. 



I — ii 



We have written 

Ai = AoUi 

where Ui is the triangular matrix with a diagonal made of 1 below: 



n + j-ii {n + j -iy. Jiei 



We now perform the transformation 



.(2) 



fc<'> 



n + j-ii 



Cf^-!^±J-1L-1-^ ifie{2,3,...,n-l}. 

The generic term of the column cj^^ , for j e {2, 3, . . . , n — 1}, is 

i — ii 1 i — ii 1 

n + j-ii {n + j-i)\ (n + j - ii)(n + j - 12) (n + j-i-1)! 

{i-ii){i-i2) 1 



0<i,i<n-l 



ifiG{0,l}, 



in + j - ii){n + j - 12) {n + j-i)\ 
This transformation supplies a matrix A2 with columns C^- , j G {0, 1, . . . , n — 1}, which writes 



(n — z)! y . 



z — zi 



1 



n + 1 — zi (n + 1 — z)! 

(Z - Zl)(z - Z2) 



(n + j - zi)(n + j - Z2) (n + j-z)! 



3e{2,3 n-1} , 



We have written 



A2 = A1U2 = A0U1U2 
where U2 is the triangular matrix with a diagonal made of 1 below: 



U2 = 6i.i - 



n + j — ii — 1 



{n + j - ii){n + j - 12) 



— 5ij_l%>2}^ 



0<i,j<n-l 



In a recursive manner, we easily see that we can construct a sequence of matrices Afe,Ufc, k G 
{1, 2, . . . , n — 1}, such that Afc = Afe_iUfc where Ufc is the triangular matrix with a diagonal made 
of 1 below: 



TT IX {n + j-ii-l)...{n + j- Zfc-i - 1) 



0<i,j<n-l 
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and 

Afe . 



I — ll 



n + l-ii {n + l-iy. J 

{i - ii) . . . (i - ik-i) 



1 



(n + fc — 1 — ii) . . . (n + fc — 1 — ik-i) {n + k — 1 — i)! J .^^ 
(i - ii) . . . (i - ife) 1 



{n + j - ii) . . . {n+ j - ik) {n + j - i)\ J iei 



We finally obtain, since all the Uk^ fc G {1, 2, . . . , n — 1}, are regular, that 

Ao = A„_iU,-ii...Ur^ =LU 



withU = U,;ii...Uj;^ and 



L = An-l 



(n — i)! y,^^ : + 1 — ii {n + 1 — i)l J 

(i-ii) ...(i- in-i) 



{2n - 1 - zi) . . . (2n - 1 - ?„_i) (2n - 1 - z)! J 

It is clear that the matrices U and L are triangular and regular, and then Aq (and A) is also 
regular. Moreover, the inverse of Aq can be computed as 

Aq-i =U-iL-i =Ui...U„_iL-i. 

The proof of Lemma [A. II is finished. □ 
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